Observation of two-wave structure in strongly nonlinear dissipative granular chains 
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In a strongly nonlinear viscous granular chain under conditions of loading that exclude stationary 
waves (e.g., impact by a single grain) we observe a pulse that consists of two interconnected but 
distinct parts. One is a leading narrow "primary pulse" with properties similar to a solitary wave in 
a "sonic vacuum." It arises from strong nonlinearity and discreteness in the absence of dissipation, 
but now decays due to viscosity. The other is a broad, much more persistent shock-like "secondary 
pulse" trailing the primary pulse and caused by viscous dissipation. The medium behind the primary 
pulse is transformed from a "sonic vacuum" to a medium with finite sound speed. When the rapidly 
decaying primary pulse dies, the secondary pulse continues to propagate in the "sonic vacuum," 
with an oscillatory front if the viscosity is relatively small, until its eventual (but very slow) disinte- 
gration. Beyond a critical viscosity there is no separation of the two pulses, and the dissipation and 
nonlinearity dominate the shock- like attenuating pulse which now exhibits a nonoscillatory front. 



PACS numbers: 46.40.Cd,43.25.+y,45.70.-n,05.65.+b 

Strongly nonlinear dissipative granular systems at con- 
ditions of loading for which no stationary solution exists 
may exhibit new unexpected behavior. In this Letter we 
present the first observation of a new phenomenon in a 
strongly nonlinear discrete viscous granular chain sub- 
jected to a very short initial impulse. The specific model 
under consideration is a chain of granules that interact 
via a power law potential, 

(1) 

V(S k>k+1 ) = o, 6>0, 

where Sk,k+i = Vk — Uk+i and yk is the displacement of 
granule k from its equilibrium position. In general the 
prefactor A is a function of the Young's modulus E, the 
Poisson ratio a, and the principal radius of curvature 
R of the grain surfaces at the point of contact. The 
exponent n depends on the topological properties of the 
contacting surfaces. For the physically important case of 
the Hertz potential, n = 5/2 (spherical granules), A — 
[£7/3(1 - <t 2 )]V2R 1]. The equation of motion for the 
fcth grain inside the chain is 

Xk = [7 (ik+i - ik) - {xk - Zfc+i)™" 1 ] 0(x k - x k +i) 

+ [7 (xk-i - x k ) + (Xk-i - Xk) 11 ' 1 ] 0(x k -i - x k ), 

(2) 



where a dot denotes a derivative with respect to t, and 7 
is the viscosity coefficient. The Heaviside function 8(x) 
ensures that the elastic and the viscous grain interac- 
tions exist only when the grains are in contact. Here 
we have introduced the rescaled variables Xk — Vk/b, 
t — (vo/b)r, 7 = ^(b/mvo) [b = (mv^/A) 1 /™} similar to 
those of Ref. Q • Note that the constant A as well as the 
mass m have been scaled out. Initially the granules are 
placed side by side, just touching each other but without 
precompression, and a velocity vq is imparted to a single 
grain (vq = 1 in the scaled problem). 

In the absence of the dissipative terms in Eq. <|2J), the 
solution to the problem is well understood as a result of 
extensive numerical and analytic studies 0,1^,0,1^- F° r 
n > 2, analytic solutions in the long- wavelength approxi- 
mation agree with numerical simulations, and agreement 
is better for smaller n. One finds that the impact by a 
single particle (equivalent to a (5-function force applied to 
this particle) quickly develops into a stationary solitary 
wave whose width/height depend on n. For the physi- 
cally interesting case of elastic spherical grains (n = 5/2) 
the solitary wave resides on about five grains. 

The different approaches to the contact dissipation can 
be found in 0, U 0> 0, El • Viscosity leads to a number 
of salient effects. In particular, dissipation based on rel- 
ative velocities ofgranis is able to dramatically change 
the pulse profile [HEi 
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FIG. 1: Breakdown of the total energy (solid line) into por- 
tions associated with the primary pulse (dotted) and the sec- 
ondary pulse (dashed) for n = 2.5 and 7 = 0.01. Inset: 
time evolution of the total energy, where two sets of three 
curves are shown, corresponding to n = 2.2 (upper set) and 
2.5 (lower set) in the Hertz potential. Each set shows results 
for 7 = 0.01 (dotted), 0.005 (dashed), and 0.001 (solid). Early 
exponential decay reflects the energy loss mainly by the pri- 
mary pulse. The energy settles to a seemingly constant value, 
which is the remaining energy stored in the secondary pulse. 
This energy dissipates on a much longer time scale than the 
initial loss. 



In our case, below a critical viscosity, a pulse simi- 
lar to a solitary wave caused by the strongly nonlinear 
dispersive forces in the discrete medium is still gener- 
ated. In subsequent discussion we shall call this pulse or 
its remnants the primary pulse. However, because the 
pulse is spatially narrow, there are high velocity gradi- 
ents that cause a relatively rapid loss of its energy. A 
quasistatic precompression appears behind the primary 
pulse because the grains now stop "too early" [Tl| as 
the pulse moves forward. This precompression is due 
entirely to dissipation and does not exist at zero viscos- 
ity. It changes the nature of the medium behind the 
primary pulse. A broad secondary pulse follows the pri- 
mary pulse. This combination is an entirely new struc- 
ture. The secondary pulse is appropriately thought of as 
a "dissipative pulse" since it only occurs in the presence 
of dissipation; it has much smaller velocity gradients than 
the primary pulse and is therefore far more persistent. It 
quickly evolves into a long-lived structure with a long 
tail of grains of uniform velocity. In the following, we 
discuss the detailed evolution of these pulses, as well as 
the behavior of the excitation above the critical viscosity. 

Starting from the impact by a single grain correspond- 
ing to an initial energy Eq = 1/2, a small amount of 
energy is lost in some back-scattering of nearby granules, 
but almost all of the energy resides in the forward travel- 
ing wave, both parts of which are formed very fast. The 
total energy of the system as a function of time is shown 
in the inset in Fig. ^ for two values of n and three values 



time 




-300 -200 



FIG. 2: Snapshots of the velocities of the grains at early 
times ranging from 40 to 95 in steps of 5 in nondimensional 
units. The abscissa is the moving variable £(t) = k — f c(t)dt, 
where c(t) is the time-dependent velocity of the primary pulse, 
and k denotes the granule in the chain. Eventually the pri- 
mary pulse and the precompression behind it vanish and the 
secondary pulse continues to move at an essentially constant 
velocity in the sonic vacuum. Parameters are n = 5/2 and 
7 = 0.005. 



of 7. The separation of the energy into the "primary" and 
"secondary" portions for the case n = 2.5 and 7 = 0.01 
is shown in Fig.^ The primary pulse is a highly nonsta- 
tionary portion of the wave that maximizes the rate of 
dissipation of some of the energy. This is reflected in the 
steep exponential decay associated with this loss. The en- 
ergy decay slows down drastically as the primary pulse 
vanishes and only the more persistent secondary pulse 
remains. Note that dissipation of energy in the exponen- 
tial decay regime is faster for higher n. We explain this 
behavior by the larger velocity gradients in the primary 
pulse, whose width decreases with increasing n 0. An 
excellent numerical fit to the primary pulse decay is pro- 
vided by the expression E(t) — Eb exp(— 0.927*), while 
the rise of the secondary pulse is A[l — exp(— 0.927*)] 
where the coefficient A = 0.005 is the maximum energy 
of the secondary pulse for these particular parameters. 
As noted earlier, the secondary pulse also decays eventu- 
ally, but far more slowly than the primary pulse, so that 
on the time scale of Fig. 2] the energy in the secondary 
pulse quickly reaches its maximum value and remains es- 
sentially constant. 

The primary pulse travels along the chain with a di- 
minishing speed since its amplitude is decreased by the 
dissipation. At very early times the secondary pulse also 
exhibits a slightly diminishing amplitude and velocity, 
but increasing energy, as it quickly settles into a broad 
pulse of almost constant velocity amplitude with a uni- 
form velocity tail that stores kinetic energy. 

Figures [21 and |3 serve to detail the behavior of the sec- 
ondary pulse while the primary pulse has not yet disap- 
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FIG. 3: Velocity profile of the secondary pulse for n = 5/2 
and 7 = 0.01 at different times, which increase upwards on the 
left side of the pulse and downward on the right, consistent 
with the steepening of the compression front of the pulse with 
time. The abscissa is now the moving variable with respect to 
the essentially constant velocity of the peak of the secondary 
pulse. 



peared. At first, the secondary pulse moves more slowly 
than the primary, but this trend reverses as the primary 
pulse slows down with its faster loss of energy and the 
peak of the secondary pulse acquires an essentially con- 
stant velocity amplitude. Both figures show that the 
secondary pulse is asymmetric, generating an extremely 
persistent tail of essentially equal velocity granules be- 
hind it (not seen explicitly in the figures). The practi- 
cally constant velocity amplitude of the secondary pulse 
is clearly seen in Fig. [21 where the abscissa for all the 
times shown is scaled with the same constant velocity c. 
It is interesting that this peak speed can be associated 
with a local "speed of sound" c s , even though this sys- 
tem is initially a sonic vacuum 0, Q and the top portion 
of the pulse is only about 10 particles wide. From the 
precompression A = xu — %k+i, where k is the position 
of the grain with maximum velocity, we obtain a disper- 
sion relation for the linearized chain whose coefficient is a 
local "speed of sound" in the long wavelength approxima- 
tion, c s = y/Jn — 1) A" -2 . The value of A was obtained 
numerically. Table [I] compares the value of c s with that 
of the speed of the peak of the secondary pulse for two 
potentials and several values of the damping parameter. 
The agreement is clearly excellent. 

Next we follow the continuing history of the pulses. In 
Fig. 0] we show a typical low-7 unsealed progression with 
time. The figure exhibits all the characteristics we have 
discussed above, but it also shows three additional fea- 
tures. One is that the secondary pulse, being nonlinear, 
continues to change in shape. The pulse steepens (be- 
coming more and more asymmetric) as its peak travels 
faster (with local velocity c s ) than the local sound speed 
at the bottom right of the peak. Secondly, the figure 
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TABLE I: Comparison of peak pulse velocity and c s calcu- 
lated from the precompression for two values of the potential 
exponent n and various values of the damping. 
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FIG. 4: Snapshots of the velocity profile for small viscosity 
(7 = 0.02) at different times whose progression is easily rec- 
ognizable as both pulses move forward, the secondary pulse 
steepens, and the primary pulse disappears. The times are 
500, 900, and 1400 and n = 5/2. Inset: detailed view of the 
crest of the velocity profile at time 1400. 



shows that the primary and secondary pulses have com- 
parable amplitudes before the primary pulse dissipates. 
Figure 0] also shows that when the secondary pulse is 
steep enough, dispersion starts to prevail and the front 
displays oscillatory motion with structures that are a few 
grains wide, similar to the primary pulse (see inset of 
Fig. The secondary pulse is shock-like, with velocities 
of the grains in the pulse at least one order of magnitude 
smaller than the pulse phase speed. 

The detailed results presented to this point are asso- 
ciated with small values of 7. In this regime it has been 
reasonable to speak of two pulses as though they were 
separate entities, the primary being mainly due to nonlin- 
earity and discreteness, and the secondary one caused by 
dissipation and nonlinearity. The primary pulse causes 
the precompression that underlies the secondary pulse, 
and in this sense both together are a single entity. Nev- 
ertheless, as a convenient manner of speaking it is not 
inappropriate for these low viscosities to speak of a "sep- 
aration" of pulses. 

For small viscosities (7 < 0.03) the secondary pulse 
reaches a critical slope before catching the primary pulse, 
while the primary pulse loses almost all of its energy be- 
fore being absorbed by the secondary pulse (Fig. . We 
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FIG. 5: Snapshots of the velocity profile for intermediate 
viscosity (7 = 0.04) at different times: 140, 220, and 400, 
with n — 5/2. 



case the excitation is caused by an impact by a single 
grain, equivalent to a (5-function force excitation. In the 
presence of weak dissipation, the total pulse consists of 
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have observed in computations not presented here that 
in this small-7 regime the maximum velocity in the sec- 
ondary pulse increases with increasing viscosity because 
larger dissipation is associated with a greater precom- 
pression resulting in a secondary pulse of higher ampli- 
tude. We also point out that for very small 7 (< 0.002) 
the secondary pulse has an almost imperceptible ampli- 
tude on our numerical scale (and of course it disappears 
entirely when 7 = 0), and the primary pulse has a very 
long life. However we do not find a transition to a regime 
without a secondary pulse for any finite value of 7. The 
secondary pulse simply fades away smoothly with dimin- 
ishing 7. 

For intermediate viscosities (0.04 < 7 < 0.07) the sec- 
ondary pulse catches up with the primary pulse while 
the primary pulse still has energy comparable to the sec- 
ondary (see Fig. [SJ. As in the previous case, after the 
first pulse disappears, the secondary pulse propagates as 
a shock-like wave with an oscillatory front caused by the 
dispersion. 

For large viscosities (7 > 0.07) there is no clear dis- 
tinction between the primary and secondary pulses. 

For very large viscosities (7 > 0.1) it is no longer ap- 
propriate to think of two separate pulses (Fig. 0. From 
the beginning, there is a single shock-like structure of dis- 
sipative origin with a sharp essentially monotonic front. 

Recently Herbold and Nesterenko found an oscil- 
latory or monotonous shock-wave structure (depending 
on viscosity) in a similar chain but with different initial 
and boundary conditions in which the velocity of the first 
particle is held constant (long input pulse). 

In summary, we have observed a new two- wave struc- 
ture consisting of a dissipation-induced shock-like wave 
and a solitary-like primary wave in a strongly nonlinear 
discrete granular chain with no precomprcssion. In our 



FIG. 6: Snapshots of the velocity profile for large viscosity 
(7 = 0.1) at different times: 100, 300, and 500, with n = 5/2. 



two qualitatively different parts. A primary pulse similar 
to that which is characteristic of a nondissipative sonic 
vacuum is formed, but it is accompanied by a secondary 
pulse whose presence is entirely due to viscosity. Because 
of its high velocity gradients, the primary pulse is rapidly 
attenuated, while the broader and smoother secondary 
pulse persists for much longer times. The velocity of the 
maximum of the secondary pulse is practically constant 
during its long lifetime, and its speed is essentially iden- 
tical to the local speed of sound. There are thus three 
distinctly separate time scales in this problem: an ex- 
tremely short scale for the formation of the double-pulse 
excitation, a fairly rapid time scale of attenuation of the 
primary pulse, and a very slow time scale for the eventual 
attenuation of the secondary pulse. Below a critical vis- 
cosity the secondary pulse develops a dispersion-induced 
oscillatory front. Above the critical viscosity, it is no 
longer possible to think of the primary and secondary 
pulses as separate entities, and the resulting excitation 
presents a monotonic front. Similar behavior can be ex- 
pected in other strongly nonlinear discrete media under 
short pulse excitation. 
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